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1 Introduction 



In studies of the behavior of a matter field on a black hole background, e.g. the possible 
evolution of the black hole by exchanging energy, charge and angular momentum with the 
field; one would like to know if the variables can be separated in the relevant field equation. 
Obviously, the study becomes much easier, if the system of partial differential equations can 
be reduced to a system of ordinary differential equations. This paper concentrates on spin-i 
(Dirac) field. 

The "no-hair" theorem in General Relativity asserts that the metrics of stationary black 
holes can be described uniquely by three parameters : Mass M, charge Qe (assuming the 
absence of magnetic charges), and angular momentum per unit mass, a. Therefore, when 
one writes a matter field equation on a black hole background, these parameters become 
parameters of the field equation. 

The Dirac equation and its separability properties on black hole metrics have been in- 
vestigated in increasingly complicated contexts. The early pioneering work was done by 
Brill and Wheeler who separated the neutrino field equation on the Schwarzschild met- 
ric. Teukolsky p| separated a master equation for massles spin 0, 1 and 2 and noted that 
the separated equations for the massless spin-i have the same form, and therefore could be 
incorporated into the master equation. Unruh [Q showed the spin-| result independently. 
The unexpected result that the massive Dirac Equation was also separable came from Chan- 
drasekhar The Chandrasekhar result was generalized to the Kerr-Newman {Q,a,M) 
case by Page p and Toop p. 

However, black holes could also have magnetic charge, if such existed. Such a black 
hole would acquire an additional label Qm for the magnetic charge. The interest in this 
possibility has grown since magnetic monopoles have been found to be required in various 
extensions of the standard model of particle physics. Dudley and Finley [§, carried out the 
separation of variables for all real, massless, single s])m field equations with s = 0, |, 1, 2 on 
the seven-parameter class of Petrov type-D solutions of Einstein-Maxwell equations found 
by Plebanski and Demianski 0|. The seven parameters include M, a, Qe and Qm (The 
coordinates are not of Boyer-Lindquist type). Since the fields are real, these results apply 
only to neutral particles. In this paper, we use the tetrad formalism to show that the Dirac 
equation for massive, charged fermions remains separable when magnetic charge is added to 
the black hole, i.e. in the dyonic Kerr-Newman case; and we present the separated radial 
and angular equations. We anticipate using this separability in a thought experiment to test 
the cosmic censorship conjecture by considering a massive, charged Dirac field on a dyonic 
black hole, i.e. the Dirac-field analogue of work done in 
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2 The Tetrad Formalism 



The equivalence principle tells us that the laws of Physics, as written in Minkowski space, 
are also valid in a freely falling frame on a curved space, i.e. a local Lorentz coordinate system 
(LLCS). This principle can be utilized to construct the tetrad formalism^, a prescription for 
writing the laws of Physics on curved spaces. We erect at each point in spacetime, four 
vectors the "basis vectors of a LLCS". Here, the latin indices from the early part of the 
alphabet label the coordinates in the LLCS, and the greek indices, the spacetime coordinates. 
The set of four Ki's, the tetrad or vierbein, is usually chosen to be orthonormal: 

9,^V/V,'' = Vat ■ (1) 

The LLCS indices are raised and lowered with the Minkowski tensor rjat (we use rjab = diag(-l, 1,1,1)) 
and spacetime indices with the metric (7^,^; the tetrad satisfies 

VatV^.V, = 9,. (2) 

We can write down components of physical fields in the LLCS: 

A^t: = v^^v'^v/v,^ ■ ■ ■ A':,;:: (3) 

These objects A'^^,,] are scalars with respect to spacetime operations, and tensors in the 
LLCS, the number of the covariant and contravariant indices being the same as for the 
original object A'^'^Z 

In terms of these components, we want to write down laws of Physics in the LLCS, 
in their Minkowskian form. We also want the laws to be invariant under local Lorentz 
transformations, i.e. different tetrad choices. This leads us to define the coordinate-scalar, 
Lorentz- vector 'derivative' 

Va = V/d, + \v/V,''V,,,,T}"^ (4) 

where S"^ are generators of the Lorentz group in the representation associated with the field 
Va is acting on; and the semicolon denotes the (metric-) covariant derivative^. 

Now the prescription is: We take the Minkowskian form of the relevant equation, replace 
each term by its LLCS equivalent via eqn. (|^) and each derivative by its LLCS equivalent 
via eqn. (HI). 



^See, e.g. Weinberg ||T§, Section 12.5, Birrell & Davies |Tl|, Section 3.8, Chandrasekhar |T2|, §7. 
•^Weinberg IM and Birell & Davies |0| give -incorrectly- an ordinary derivative here. 



2 



3 The Dirac Equation on the dyonic black hole metric 

The Dirac Equation in flat space with a background electromagnetic field is 

z-f^'ida + ieAa)^ = fl^ (5) 

where \Ef is a four-component spinor, /i the mass and e the charge of the field quantum, A 
the vector-potential of the electromagnetic field and the constant 7-matrices satisfy 

{7^7^} = -2r/'^^ (6) 

We will be using 7-matrices of the form 

^-(:. 0). ^-(!-r) 

where the a-matrices are the identity matrix and the Pauli matrice s: 



1 \ 1 / 1 \ 2 / -M 3 1 
"^=i 1 ' " = 10' " =U ' " = -1 



.0 



For a spin- 1 field, the S-matrices are given by 

S^'^ = -^[?,7l. (9) 

To get the Dirac Equation in curved space, we follow the prescription and make the 
replacements (0,|D: 

^bVV^'d,^, + l^bVV/V,^'V,.,, [J:'V^, - eWV,>^A,)^, =^^^, (10) 

where we have written out the spinor indices j, k, I explicitly. 
The metric of the dyonic black hole isQ 

ds^ = g^^dx^'dx" = df + ^ sin^^ rf^^ 

^2 

ubin u uiu(p -\ — —ui' -|- /"^'"^^ 

P 

where 

p2 = r^ + a^cos^^ (lib) 
A = r^-2Mr + a^ + Q^^ + Q^^ (11c) 



2^-^4^^asin2^ dtdcP + ^dr' + p^d9^ (11a) 



*For derivations in Boyer-Lindquist coordinates, see H and ||13|. 
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with the vector potential 



^ r ^ acosO 
A = -Qe^ + Qm — T- 12a 

= Ae = (12b) 

A4, = Qe 5 + g„±l-COS0 —]. (12c) 

P"^ P"^ 

and we are using Boyer-Lindquist coordinates. 

Of course, the vector potential is unique only up to a gauge transformation, and the 
magnetic part of A^ contains a string singularity. The two signs in that term correspond 
to the two gauges we will be using. The upper-sign term puts the string along the negative 
z-axis {6 = tt) and will be used when < < 7r/2, the lower-sign term puts it along the 
positive z-axis {6 = 0) and will be used when n /2 < 6 < n. Therefore, the wavefunction is 
also gauge transformed across the equator, and picks up a factor of e^*^*^™'^ passing from 
north to south. This matching of boundary conditions ensures that the problem can be 
expressed meaningfully without strings of diverging vector potential. Such a wavefunction 
is called a section lO . 



We will derive separated equations first using a simple tetrad that we call 'canonical'. 
For purposes of comparison at the appropriate limit with Chandrasekhar Q], we will repeat 
the calculation using a more complicated tetrad, calling it 'Kinnersley-type'. The labels (C) 
and (K) will be used for variables whenever it is not obvious which case they belong to or 
for emphasis. 

3.1 The Canonical Tetrad 

We use the generalization of the simple 'canonical' tetrad emphasized by Carter, relabel- 
ing the vectors to be able to compare with Carter & McLenaghan |jl9|] in the limit Qm = 0: 





-{dt 




P 




sin^^ 




P 




pde 




P-dr 








- a sin 6 dcj)) (13a) 

-a dt + (r^ + a^) dcf))] (13b) 

(13c) 
(13d) 

where 5 = VA and p was defined above via its square. 

Making all the necessary substitutions into eq. (p!0D , we get the four coupled components 
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of the Dirac equation^ on the dyon black hole metric: 

5 „ a „ 6 



po p po 

a sinO _ i „ 

+ dt + -de 

P P 



2pp* 



(r-M) , , -Qer±aQr 



pS 



1 ^ a sinO ^ i . ^ ±1 — cos^ 

psinO 2pp* 2p ™ psinO 



-ip'^P (14a) 



a sinO „ i „ 1 „ a sin6' i . ±1 — cos^ 

dt - -de —d^ + — — ^ - — cot^ - teQm r^— 

p p psmu 2pp* 2p psm& 



+ 



(r2 + a^' 



p5 



dt + -dr + 



P 



p5 



6 (r-M) -QeT ± aQr 



2pp* 



2pS 



dt + -dr + — ( 

po p po 



5 (r-M) -Q^r ± aQr 
+ — + ^ . ^ +ie- 



2pp 2p6 



p6 



pS 

(c) 



+ 



asinO i 1 asinO i n ■ ^ ±1 — cos6' 

dt de H —d^ — — - — cote* + leQm r-^— 

p p psmtl 2pp 2p psm& 



asinO ^ 1 „ asmO i n ■ ^ ±1 — cos6' 

dt + -de H —d^ + — — + — cota* + leQm 

p p psmU 2pp 2p psmtl 



(C) 



+ 



(r^ + 



_ ^ _J__ (r-M) ^ .^ -Qer±aQr 

pS p p6 2pp 2p6 p6 



^1 



i/i^P (14b) 



(14c) 



-ip'^P (14d) 



where 



p = r + ia cos9, p* = r — ia cos9, pp* = p^ 



3.2 The 'Kinnersley-type' Tetrad 

For sake of comparison with Chandrasekhar @], in the limit Qe = 0, Qm = 0, we write 
down the Dirac Equation using the generalization of the Kinnersley-type tetrad correspond- 
ing to the null tetrad used by him: 



Y(k)o 

Y(K)3 



(1 + -^) dt + (-- ^ dr 
^/2^^ 2p^' 4 



'1 



A 



cos6 sinO 
? 



dt + r dr + 



a(a2 + r2)cos^sin^ ,^ 
5 d(p 

p2 



arsin6' , ^ ,^ r(a^ + r^) sin6' 
— dt + a cosO de - — 



P 



P 



1 



A 



lasin 6 



(15a) 
(15b) 
(15c) 
(15d) 



^The calculations for going from eqn. ( [lO| ) to eqns. (14a- 14d) or eqns. ( [I6a| - |l6d ) have been performed 
using the Mathematica symbolic mathematics software, and the 'tetrad' package ]15[. 
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Once again, making necessary get four coupled components of the Dirac equation on the 
dyon black hole metric; for the Kinnersley-type tetrad: 



+ 



-dr H T^-ra 



r — M) ^ ie f ±aQm — Qef' 



ia sm9 1 i 
— —dt - -de - — 
p p p smt^ 



ia^voB 1 ^\ ^ /±1 — cos^^' 

— + — cot ^ ) + tQr: ' 



p 



2-p 



ia sine ^ _ ^ i q Lcot6'-eQ | ±1 - cos^ 

p* * p* ^ p* sin6' 2p* ™ y p* sin6' 



y p sin6' 



ip'^f^ (16a) 



+ 



A A p* V A 



-zp^S^^ (16b) 



A A p V A 



(i^) 



+ 



msin6'„ 1^ ^ o 1 r, ^ ( ±1 — cosO^ 
— - — dt + -de + —^04, + —coiO - eQm — _ . ^ 
p p psmt/ 2p Y psmt7 



[K) 



ip^f^ (16c) 



ia sm9 „ 1 ^ i 
— —dt + -de - --— 
p* p* p* smu 



iasmO 1 ^\ ^ /±1 — cos6'\ 
— + — cot ^ ) + eQ,-- I I 



P 



2p* 



\ p* sinO J 



(K) 



+ 



(r^ + a^] 
V2p' 



A a r - M ie I ±aQm - Qer 
dt T^^dr + ^ ^ d^ + 



V2 



p- 



V2 



p- 



y2p2 V2 



p- 



ip^f^ (16d) 



4 The Separation 



To separate eqns. (|14a| - [L4d|) or ( |16a| - |16dD , we assign the standard time- and azimuthal 
dependence to \E', with the above-mentioned gauge transformation across the equiatorial 
plane: 



= e 



-iujt i{m^eC, 



(17) 



4.1 The Canonical Tetrad 



Substituting eqn. (^) into ( |14a| - p.4d|) gives, after some algebra 



6 ,^ 1 , , i . ^ ia sind . , 
p 2p* p 2p* 

iasm6 , , 5 1 



P 



2p* -)^^ + p(^- + 2^)^^ 



ia sinO 



2p 



5 1 X , i , ^ 

p 2p p 

i , „ ia sm.6 . , 5 , 1 , , 

P 2p p 2p 



-ip^o 
—ipipi 

-ipip2 
-iptps 



(18a) 
(18b) 
(18c) 
(18d) 
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where T>^£^ and C±^^ are purely radial and purely angular, respectively: 



(r^ + 0^)00 + eQeT — ma (r — M) 
Or + i : h 



A 



2A 



V 



.{r"^ + 0^)00 + eQ(.r — ma {r — M) 
~ ' A + 2A 



(19a) 
(19b) 



de — auj sinO + 



m 



sm9 



{- + eQ^)cot9 



TTl 1 

de + auj smd + ( eQm) cot 9 

sm.0 2 



(20a) 
(20b) 



Defining = /o/Vp, ^1 = fi/y/p, = and ^3 = /s/V/^; multiplying eqns. ( [l8aD 

and (|18bD by pi/p*' eqns. (|18c| ) and (|18d| ) by Pa/p, we get 



-5 V^h+iC+f, 



-ifi{r — ia cos9)fo 
-ifi{r — ia cos6)fi 
-i^{r + ia cos9)f2 
-ifi{r + ia cosO)/^ 



(21a) 
(21b) 
(21c) 
(21d) 



The structure of eqns.( pTa| - |21d| ) suggests that they can be separated by the substitutions 

/o = R-S+ ; /i = R+S- ^22-) 
/2 = R+S+ ; /s = R-S- 

which gives 



6 T>^R^ — ifirR^ = XqR_ 
6 V_R_ + ipri?+ = Aii?+ 
5V^R^+ ifirR+ = X2R+ 
6 V+R+ - i^rR_ = X-^R- 

The consistency of these equations requires 



iC^S^ + fiacos9S+ = XqS+ 
iC^S^ — fiacosOS^ = XiS^ 
iC-S^ + fiacosOS-^- = A25'+ 
iC^S+ — ^acos9S_ = X-^S^ 



-^0 — -^3 — -^1 — -^2 — 



and we are left with 



,(C) d(C) 



ifirR 



S V^y'R^L' + ifirR 



(23a) 
(23b) 
(23c) 
(23d) 



(24) 



(25a) 
(25b) 
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and 

Cf'sf^ + ifia cose S^^^ = -iXS^^^ 

Or, by combining pairs, we can get decoupled 2'^'^ order equations 

[6 - --^^Pf - (A^ + /iV2)]i?f = 

[6 Vf^6 P(^) + - (A^ + A')]^^""^ = 
(A — ifir) 

[^(C)^(C) ^iasin^^(c) ^ _ 2^2 ,^,2Q^^sf^ ^ q 

^ + A + /iasin6' + ^ + 

X — fia cose* 

4.2 The Kinnersley-type Tetrad 

The same time- and azimuthal dependence ([T7|), when substituted into eqns. 
gives 

A ^ , 1 , - ia sm9 . 

+ H =;^)V'3 = Wo 



V2{v. + 4)^3 - 4-^-^2 



72 = -^PV'I 

p" ■ p" 



-«P^2 

P' P 

A 1 iasin^ . 

^+^1 - ^(^- —m = «P^3 



v^p2 P*' p 

the radial and angular operators and £fif ■* being 



= dr -^[{r'' + 0^)10- ma + eQer] 



C^+ ^ = dg + au sinO + + eQm) cot 6 

smd 2 

C[^^ = de — auj sinO H + ( eQm) cot 6 

sm9 2 



,2 _ ^„ '' '29 bp by p*, eq. ( p9cD by p; defining 
-^0 = pV^o, = ^Ji, F2 = ^2, and F3 = p*V^3, we get 



Multiplying eqns. ( |29a| ) and ( |29d| ) by = pp", eq. 



\/2P_F3 - /:_F2 



Again, the eqns.( ^2a| - |32d|) can be separated by 



F2 



ip{r — ia cos9)Fq 

—ip{r — ia cos9)Fi 
—ip{r + ia cos9)F2 

ip{r + ia cos6)F3 



Fi — R^S— 
F3 = R_S^ 



which gives 



as before, 
therefore 



V2 



iprR 



^vi^)R(_^) + ,^rRi^) = x'R 



(32a) 

(32b) 
(32c) 

(32d) 



(33) 



^V+R+ 


— iprR- = 


XqR- 






— pa cos9S+ 


— — Ao5'+ 


(34a) 


y/2V.R^ 


+ iprR+ = 


A1-R+ 


; 


-S+ 


— fia cos9S- 


= Ai5_ 


(34b) 


V2V^R^ 


+ iprR+ = 


A2-R+ 


; 




— pa cos9S-^- 


= — A25'+ 


(34c) 


-^V+R+ 
\/2 


— iprR_ = 


As^?- 


; ^- 


-S+ 


— pa cos9S^ 


= A35_ 


(34d) 




Ao 


= A3 = 


Ai = 


A2 


= X' 




(35) 



(36a) 
(36b) 



ci^^S^^^ - pa cos^^f ) = -A'5f ) (37a) 
ci^^sP - pa cos9Si^^ = X'si^'^ (37b) 

The corresponding decoupled 2'^'^ order equations are 

lAvi^)jyW + //^^ p(^) _ + ^2r')]R(^^ = (38a) 
(A' — ipr) 

[p(^)Apf ) - Apf) _ (A'2 + pV2)]i?f ) = (38b) 

(A' + ipr) 
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^ + X' + fia cose V A- yj + 



(i^)^(i^) 



X' — na cos9 








(39a) 
(39b) 



5 Comparison with previous results 



Both papers that we compare our results against use the Newman-Penrose formahsm. For 
the simpler 'canonical tetrad', our results (eqns. |25al - p8b|) agree with Carter & McLenaghan 
1| in the limit = 0, with V^i, Cf^ ^ X^i, Sf '' 



Yi. 



S. 



-r 1, A 



-\/2 A, yU v^/i. For the Kinnersley-type tetrad, the equations ( p6a 



|39b| ) agree with Chandrasekhar []T2| in the limit e = 



Qm 0, 



with T> 
-R 



X' 



2 

-A. 



0, 

2 ^ ^ ^ ^ 

A difference with the Qm = case is that m can now also take on half-integer values: To 
make the wavefunctions ([TtD single- valued, both m + eQm and m — eQm have to be integers, 
which makes m and eQm integers or half-integers, the latter requirement being, of course, 
the well-known Dirac quantization condition. Since the explicit forms of equations ( P7a| - |28b]) 
and ( p8a| - |39b]) are not very illuminating, we leave them in their present form. Unfortunately, 
the angular equations ( |284|28b| , |394|39b| ), unlike in the simpler scalar case |]T6[, are not of 
the Sturm-Liouville form, therefore we cannot make a statement about the completeness of 
the solutions. 
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the 



The Chandrasekhar result has been also generalized to other cases 
mathematical structure of the problem has been investigated , |T^ , , |^ , , |^ , ,P5| 



Other approaches to the problem of separation of 



and the solutions studied 
variables in the Dirac Equation on curved spaces include the Stackel Space method [0 and 
the "algebraic" method pl| 
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